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Abstract 

The cascade model generates random food webs. The contin¬ 
uum cascade model is a Poisson approximation of the cascade model. 
We have a simple recursion to obtain probability distribution of the 
longest chain length (the height) generated by the continuum cascade 
model. Assuming the traveling wave solution, the velocity selection 
principle for the Fisher-KPP equation works to the recursion, 

Here we formulate the height of continuum cascade model as the 
first passage time of the left most particle of a branching Poisson 
point process. The asymptotic probability distribution of the height 
is obtained mathematically by a straightforward application of the 
Aidekon theorem for the left most particle of branching Poisson point 
process. 


1 Introduction 

Random graphs are modeling numerous natural phenomena, polymerization, 
spread of infectious diseases, transportation systems, electrical distribution 
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systems, the Internet, the world-wide web, social networks, etc. ( see [23] and 
references therein). A random graph is a set of vertices that are connected 
by random links. In random graph modeling, links are usually treated as 
undirected. However, directionality plays a prominent role, in modeling of 
the web growth, modeling of food webs, etc. The cascade model ( [11, 10]) is 
introduced to analyze the ecological data for community food webs and is a 
standard model in the subject. The model generates a food web at random. 
In it the species are labeled from 1 to n, and arcs are given at random between 
pairs of the species. For an arc with endpoints i and j (i < j), the species i 
is eaten by a species labeled j. The longest chain length (the height) is 
compared with the ecological data [11]. The cascade model provides a natural 
mechanism for generating directed random graphs and the same model has 
been suggested in other contexts, e.g. as a model of parallel computation 
[17] where the presence of the directed link (i, j) with i < j indicates that 
task i must be performed before task j. For a parallel computation in which 
each task takes a unit of time, the processing time will be equal to Mn -l- 1 
for the length of the longest path M„. 

The continuum cascade model is a Poisson approximation of the cascade 
model [23]. We have a simple recursion to obtain probability distribution 
of the longest chain length (the height) generated by the continuum cascade 
model. We applied the velocity selection principle [24] for the Fisher-KPP 
equation to it assuming the traveling wave solution [23]. The Fisher-KPP 
equation [15, 24], which is for the spacial spread of an advantageous gene, 
has a traveling wave solution with the velocity determined by the velocity 
selection principle. The branching Brownian motion gives a solution to the 
Fisher-KPP equation and gives the velocity of the wave front with the cor¬ 
rection term [7, 26, 27, 34]. 

Here we obtaine the recursion for the probability distribution of the height 
of continuum cascade model [23] from the position of the first passage time 
of the branching Poisson point process. We have the asymptotic probability 
distribution of the height mathematically by a straightforward application 
of the theorem [4, 5] on the left most particle of the branching random walk 
[2, 20]. The Aidekon theorem is the branching random walk version of the 
Lalley and Selke theorem for branching Brownian motion [26] and described 
in terms of a functional of the limit of the derivative martingale associated 
to the branching random walk. 

Our recursion (2), given in [23], seems to be very typical like the Fisher- 
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KPP equation. Actually it is pointed out that our continuum cascade model 
is also studied under the name of Poisson weighted inhnite tree (PWIT, [1, 3, 
16]). The closure of vertex 1 in the cascade model converges in distribution to 
the PWIT as n tends to inhnity [16]. The expected position of the wave front 
and the hnite width of the wave front for our continuum cascade is obtained 
for the the PWIT mathematically [1], independently from the mathematical 
study [4] and independently from our intuitive physical study [23]. 

Extending the argument on 1-dimensional random sequential packing 
([30]), we have the random sequential bisection of intervals [33] (continuum 
binary search tree), which has an analogous asymptotic behavior to the bi¬ 
nary search tree of n keys for the sorting algorithms [13, 31]. The nonlinear 
recursion for the probability distribution of the minimum of gaps generated 
by 1-dimensional random sequential packing [21, 22] is developed to the non¬ 
linear recursion for the probability distribution [31, 33] of the height of binary 
search tree. Assuming the traveling wave solution for the Hattori and Ochiai 
conjecture [18, 19, 25], the correction term in the asymptotic expected height 
of continuum binary search tree is obtained in [25]. The correction term 
is also obtained mathematically [14, 29, 32] without using the assumption. 
Like the case of binary search tree, making continuum model gives a natural 
mathematical arguments for the original cascade model. 

2 A recursion for the continuum cascade model 

In the cascade model, the random directed graph has vertex set {l,...,n} 
in which the directed edges (i, j) occur independently with probability c for 
i < j and probability zero for i > j. Let denote the length of the 
longest path starting from vertex 1. We apply the Poisson approximation to 
the binomial distribution of the number of directed edges at each vertex and 
consider the continuum cascade model. We study the probability distribution 
of Ln as n tends to infinity. [10] [28]. 

At step 1 we generate points by the Poisson distribution Pr{N = 
k) = on [0,a;]. Each point is mutually independently distributed uni¬ 

formly at random on [0, x\. At step j(> 1), for each generated point at x — y, 
generated at the step j — 1, generate Ny points by the Poisson distribution 

k 

Pr{N = k) = uniformly at random on the interval [x — y,x\, indepen¬ 

dently from the points on other intervals at step j and independently from 
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the points of previously generated intervals. We call the terminal interval 
any interval which did not generate any point. The terminal interval, which 
appeared at a step, remains a terminal interval after the step. We continue 
the steps as long as we have at least one interval which is not a terminal 
interval. For each terminal interval generated by the above procedure, we 
count the number of steps to get the interval. Let us call the maximum of the 
numbers, H{x), the height of the tree generated by the continuum cascade 
model on [0, x]. 

When k points, x — yi,x — y 2 ,...,x — yk are generated at step 1, the proba¬ 
bility, that the height is not larger than n—1, is Pn-i{yi)Pn-i{y 2 ) ■ ■ ■ Pn-iijjk)- 
Since each yi is distributed uniformly at random on [0, x] and k is distributed 
by the Poisson distribution, we have the following recursion [23] for the prob¬ 
ability Pn{,x) = P{H{x) < n). 

For n = 0, 


Pn{x) 


= e 


( 1 ) 


while for n > 1, 


Pn{x) 


1 pX pX 

e-" + E7Te""-/ ■■■/ Pn-i{yi)---Pn-i{yk)dy,---dyk 
^ k\ x^ JO Jo 


oo 1 / 

^ Pn-i{y)dy"^ 


exp 

—X -|- 

[ Pn-i{y)dy) . 

Jo J 

(2) 


3 Traveling wave 

We apply the Aidekon theorem [4, 5] for branching random walk to show 
mathematically the following observations [23] for equation (2) in later sec¬ 
tions. It is convenient to think about x and n as space and time coordinates, 
so that the front of the traveling wave was advancing 

1. Numerical traveling wave solution. Numerically, the probability 
distribution Pn{x) has a traveling wave shape with the width of the 
front remaining hnite as shown in Fig. 1, which means the hnite width 
of the height distribution. 
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Figure 1: Traveling waves, the distribution Pn{x) versus x obtained by iter¬ 
ating the recursion. Pn{x) is shown for n = 20,40, 60, 80,100 (left to right). 
Iterations were performed for a discrete approximation (29) given in section 

6 . 

2. Velocity selection. Assuming an approximation by a traveling wave 
form for larger n 3> 1, 

Pn{x) ^Ii{x - Xf), (3) 

with the front position Xf growing linearly with ‘velocity’ equal to e~^ 
[23]: ^ 

Xf ~ vn, V = -. (4) 

The velocity selection principle [24] gives v = e~^ [23] by using an 
analogous argument to the case of binary search tree [25]. We see that 
the wave front Xf should advance asymptotically by a constant velocity 
V = e~^, from the probabilistic argument for the cascade model [10, 28]. 

3. Logarithmic correction. An analogy [7, 34] from the Fisher-KPP 
equation gives a logarithmic correction to the front position as 

T) S 

Xf = —I-Inn-I-O(l). (5) 

e 2e 

It is convenient to think about x and n as space and time coordinates, 
so that the front of the traveling wave was advancing. 
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4 Branching random walk for the solution of 
the recursion 

We consider the asymptotic behavior of branching random walk [4, 5, 2, 8, 12] 
[20]. We follow the notation and argnment by Aidekon [4], The process starts 
with one particle located at 0. At time 1, the particle dies and gives birth to 
a point process £. Then, at each time n E N, the particles of generation n 
die and give birth to independent copies of the point process £, translated 
to their position. If T is the genealogical tree of the process, we see that T 
is a Galton-Watson tree, and we denote by |a;| the generation of the vertex 
X E T (the ancestor is the only particle at generation 0). For each a; G T, 
we denote by V{x) E R its position on the real line. With this notation, 
(G(a;); |a;| = 1) is distribnted as C. The collection of positions {y{x)-^x E T) 
dehnes onr branching random walk. 

We assnme that we are in the bonndary case [20] 

B[i: ii>i. 

hl=i 

E\Y. e-'''0 = l. 

hl=i 

E[Y. V{x)e-^^^'^] = 0 

|x| = l 

Every branching random walk satisfying mild assnmptions can be rednced 
to this case by some renormalization. Notice that we may have 

^ 1 = oo (9) 

hl=i 

with positive probability [5]. We are interested in the minimnm at time n 

:= min{E(a;); |a;| = n}, (10) 

where min |0| = oc. Writing for y E i?U{±cxD}, y+ := max(|/, 0), we introdnce 
the random variable 

X ;= ^ (11) 

hl=i 

X ;= ^ V{x)+e-^^^\ (12) 

hl=i 


( 6 ) 

(7) 

( 8 ) 


6 



We assume that the distribution of C is non-lattice, we have 


E[Y,V{xfe-^^^^]<oo (13) 

|x|=l 

E[X{\n+Xf] < oo, E[X{\n+X)] < oo (14) 

To state the result, we introduce the derivative martingale, dehned for any 
n > 0 by 


Dn := E (15) 

\x\=n 

From [5, 6] ( Proposition A.3 in the Appendix [5]), we know that the mar¬ 
tingale converges almost surely to some limit Doo, which is strictly positive 
on the set of non-extinction of T. Notice that under conditions (6), (7), (8), 
the tree T has a positive probability to survive. 

There exists a constant C* > 0 such that for any real x, 

lim P{Mn > +x) = (16) 

(Theorem 1 in [4, 5], see [9] for an elementary approach). 


5 Probability on the longest chain length 

The process starts with one particle located at 0. At time 1, the particle 
dies and gives birth to the point process £, with intensity 1 on [0, cxd). Then, 
at each time n E N, the particles of generation n die and give birth to 
independent copies of the point process C, translated to their position. At 
each time , we kill all particles to the right of x. Denote position of left-most 
particle in this (extended) tree at n-th generation by i7„. Since 

P{H{x) <n-l) = P{Hn>x), (17) 

we see 

Pn-i{x) = P(i7n > x). (18) 
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To normalize for equation (8) we replace the original Poisson Point Process of 
intensity 1 on [0, cx)) by the Poisson Point Process of intensity - on [—1, cx)), 
as £ in section 4. Then the conditions (6), (7), (8) hold. We see the inequality 
(6), since expected number of children here is infinite. For the identity (7) 
we have 



and for the identity (8) we have 


_dy 


1-1 


ye 


= 0 . 


(19) 


( 20 ) 


The distribution of this Poisson point process C is non-lattice and the mo¬ 
ment conditions (13) and (14) hold by exponential decay (7) for V{x). We 
have 


E[J2 

\x\=l 

/ OO 

x^e~^dx < OO. (21) 

The total number of children is assumed to be hnite almost surely in [4]. 
However the argument [4] is applied to the above extension to [— 1,cxd), as 
shown in [5]. 

The position of the left-most particle at generation n is given by using 
M„ for 


M„ = eHn — n.. (22) 

Considering equation (22), 

Mn>z + ^\nn, (23) 


means 


z + n + ^\nn 

Hn 


e 


(24) 




Figure 2: Pn-i{^ Inn + ^) for n = 1, 2,100 


Hence from equation (18), 


^ ^ 3 Inn, 

P[Mn > z --—) 


, z + n + ^ In n , 

= P{Hn > -^- 

e 

+ n + I Inn, 

= Pn-i{ -^- . 


(25) 

(26) 


Put z/e = X, then from equation (16) ( Aidekon [4] [5]), for the solution Pn-i 
to equation (2) we have 


71 / 3 

lim Pn-iix H-h — Inn) = i?[exp(-C'*e®"'Zloo)], (27) 

n—>-oo ^ g 26 L i \ /j' \ / 

which gives the asymptotic probability on the longest chain length (on the 
position of wave front). 

For a; = 0 of equation (27), we have 


71 / 3 

lim Pn-i{- + — Inn) = E[exp{-C*D^)], (28) 

n^oo e 2e L r-v /j7 

which should be less than 1 and larger than 0, since is mathematically 
shown to be strictly positive [4, 5, 6]. 
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6 Numerical observation 


Putting xA for x, and giving the discrete initial value for equation (1), we 
consider a recursion as a discretization of equation (2) for 1 < n, 


fn{x) = exp[-xA + Y, fn-i{y)A]. 

y=l 


(29) 


The numerical value fn{x) for Pn{x) in Fig. 1 and Fig. 2 are obtained 
from equation (29) for A = 0.01 by using the software Mathematica. The 
numerical values 


-)) 

A 2e e 


(30) 


for P„_i(A Inn + ^)), are shown by the lower curve in Fig. 2. Putting eja 
instead of exponential e, the numerical values 


a 


fn-ii-ri— Inn 

•' ^AAe 


n 


)) 


(31) 


for A = 0.01 and a = 0.9855 are shown by the upper curve in Fig. 2, which 
seems to approach quickly to a constant. We carried out calculations and see 
for example A = 0.001 and a = 0.9977 the value of (31) quickly approaches 
to a constant. Our numerical calculations seem to suggest a —)■ 1 as A —)■ 0, 
which supports equation (28) for the wavefront numerically. 
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